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ABSTRACT 



Context. The average Doppler shift shown by spectral lines formed from the chromosphere to the corona reveals important information 
on the mass and energy balance of the solar atmosphere, providing an important observational constraint to any models of the solar 
corona. Previous spectroscopic observations of vacuum ultra-violet (VUV) lines have revealed a persistent average wavelength shift 
of lines formed at temperatures up to 1 MK. At higher temperatures, the behaviour is still essentially unknown. 
Aims. Here we analyse combined SUMER (Solar Ultraviolet Measurements of Emitted Radiation)/SoHO (Solar and Heliospheric 
Observatory) and EIS (EUV Imaging Spectrometer)/Hinode observations of the quiet Sun around disk centre to determine, for the 
first time, the average Doppler shift of several spectral lines formed between 1 and 2 MK, where the largest part of the quiet coronal 
emission is formed. 

Methods. The measurements are based on a novel technique applied to EIS spectra to measure the difference in Doppler shift between 
lines formed at different temperatures. Simultaneous wavelength-calibrated SUMER spectra allow establishing the absolute value at 
the reference temperature of T a; 1 MK. 

Results. The average line shifts at 1 MK < T < 1.8 MK are modestly, but clearly bluer than those observed at 1 MK. By accepting 
an average blue shift of about (-1.8 ± 0.6) km s"' at 1 MK (as provided by SUMER measurements), this translates into a maximum 
Doppler shift of (-4.4 ±2.2) km s~' around 1.8 MK. The measured value appears to decrease to about (-1.3 + 2.6) km s"' at the Fe xv 
formation temperature of 2.1 MK. 

Conclusions. The measured average Doppler shift between 0.01 and 2.1 MK, for which we provide a parametrisation, appears to be 
qualitatively and roughly quantitatively consistent with what foreseen by 3-D coronal models where heating is produced by dissipation 
of currents induced by photospheric motions and by reconnection with emerging magnetic flux. 

Key words. Sun: corona - Sun: transition region - Sun: UV radiation 



1. Introduction 

Ultraviolet observations of the Sun near disk centre show a 
systematic net red shift of emission lines from the transi- 
tion region (TR). Several vacuum ultra-violet (VUV: 16 to 
200 nm) spectrometers with different spatial resolutions such 
as HRTS and those flown on Skylab, OSO 8, and SMM have 
revea l ed this phenomenori since the 1970s (e.g. , Doschek et al. 
[19761 iHassler et all [19911 iBrekkd [T99l [Achouret all Il995l) . 
Systematic red shifts h ave also been o b served in spectra 
of late type stars (e.g., lAvres et al.1 119881: IWood et all [l996t 



iPagano et al.l |2004|) . indicating that this phenomenon is not 
unique to the Sun. 

In these earlier investigations (before the launch of SoHO 
in 1995), the magnitude of the red shift was found to increase 
with temperature, reaching a maximum at T = 0.1 MK and then 
decreasing toward higher temperatures. However, these observa- 
tions were, with few exceptions, restricted to temperatures below 
0.2 MK. 

Observations from SOHO extended the observable temper- 
ature range up to about 1 MK. In the quiet Sun near disk 
centre, Doppler shifts in the range of 10 to 16 km s ' were 
confirmed in lines formed fr om T = 0.1 MK to 0.25 MK 
(e.g.,|Bre kke et al. j l l997l: [aiae'etal][T99M iPeter & Judgdll999t 
iTeriaca et al.a999alK 



At higher temperatures (0.6 to 1 MK), the int erpretation of 
the me asur ements prov ed to be more controversial. Brek ke et al.l 
( Il997h and lChae efaD ( 1998) reported red shifted (downward) 
average Doppler velocities of 6.0 and 3.8 km s"' for Mg x A625 
(r = 1.0 MK) and of 5.0 and 5.3 km s"' for Ne viii ^770 
(T = 0.63 MK). These Doppler velocities were derived by us- 
ing the refere nce rest wa velengths of 770.409 A and 624.950 A 
compiled bviKellvl ( 1987^ (based on solar observations, see dis- 
cussion in iPeter & Judgeiri999l) . affected by high uncertainties 
(high in relation to the Doppler shifts under discussion). 

Seter&Judgi (1999) and iDammasc h et al., (1999ab) took 
another approach and obtained accurate rest wavelengths for the 
above lines by assuming that the average Doppler velocity is 
around zero above the solar limb (motions along the line of sight 
cancel out on average in an optically thin plasma). With these 
new rest wavelengths, average blue s hifts of a few km s~' were 
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found for both lines near disk centre (lPeteilll999HPeter & Judgj 
1 19991: IDammasch et al.lll999a|[bl:lTeriaca et al.lll999ab~ These re- 
sults suggest that there is a transition from red to blue shifts 
above 0.5 MK. This transition is an important observational con- 
straint to distinguish between different models of the solar tran- 
sition region. 

Several hypotheses to explain the observed shifts of TR and 
coronal lines have been discussed in the literature. Since spicules 
carry an upward mass flux about 100 times larger than car- 
ried out by the global solar wind, it was suggested that the 
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TR emission is dominated by down flowing, cooUng plasma 
injected into the corona by spicules (iPneuman & Kopplll978t 
[Athay & Holzer 1982; Athay 1984). Lack of direct observational 
evidence (IWithbro e 1983) and failure of theoretical models to re- 
produce observations ( Mariska.l987J have led to a general dis- 
missal of this hypothesis. However, recent analysis of spectro- 
scopic and imaging VUV data with high sensitivity and spatial 
resolution has revived interest in the role of spicules in th e mass 
and energy balance of the solar coro n a (e.g., i De Pontie u et alj 
l2009HMcIntosh & De Pon'tieull2009allbHDe Pontieu et al.il201 ih . 

Much effort has also been put into modelling the effect of 
energy deposition (e.g., mimicking the effect of magnetic recon- 
necti on in di f ferent regions of the solar atmosphere) on line pro- 
files. IChend (Il992h showed that impulsive energy release in the 
photosphere generates a shock train that propagates upward and 
lifts and heats the chromosphere and TR, giving rise to an aver- 
age downward material velocity. However, iHansteen & Wikstoll 
(11994) demonstrated that, contrary to observations, the result- 
ing line emission is neverthele s s blue shifted. Instead of upward 
propagating waves, iHansteeril ( Il993h interpreted the observed 
red shift as caused by downward propagating acoustic waves. 
Nanoflares releasing relatively small amounts of energy near 
the top of the magnetic loops give rise to distur bances running 
downw ard along the loop legs. In a later studv IHansteen et aU 
( 119971) included the effect of reflection in the chromosphere 
of downward propagating MHD waves. This increases the red 
shifts in the TR, alt hough coronal blu e shifts become too high. 
An investigation bv lJudge et al.l (Il998h lends observational sup- 
port to the picture of downward running waves or disturbances. 

More generally, models of energy deposition at TR tem- 
peratures in 1-D loops have managed to reproduce the ob- 
served TR red shifts and low corona blu e shifts relatively well 
dTeriaca et al.lll999bl:ISpa"daro et al.ll200 6|). 

Recently, three-dimensional (3-D) comprehensive numerical 
models of the solar atmosphere (from photosphere to corona) 
have been d eveloped. In the model described by iPeter et aH 
(12004 I2006h . coronal heating is caused by the Joule dissipa- 
tion of the currents that are produced when magnetic field is 
stressed and braided by photospheric motions. This model does 
produce red shifts at TR temperatures but also in the low c orona , 
where blue shifts are instead observed. IHansteen et alj (1201 Ol) 
extended that work and added the effects of episodic injections 
of emerging magnetic flux that reconnects with the existing field. 
Rapid, episodic heating of the upper chromospheric plasma to 
coronal temperatures naturally produces downflows in TR lines 
and slight upflows in low coronal lines. As suggested by the au- 
thors, the localised heating events occurring in the chromosphere 
generate high-speed flows in a wide range of temperatures that 
may be related to the high-speed upflows that have been de- 
duced from significant blueward asymmetrie s in EIS/Hinode 
observations of many TR and c oronal lines (Haraetal. ''200Jt 
iDe Pontieu et al.1 120091: iMcIntos h & De Pontieu 2009a b) that 
have been linked to spicules (i.e., De Pon tieu et al.ll201 li) . 

With the launch of the Extreme Ultraviolet I maging 
Spectrograph (ElS. ICulhane et alJl2007l;lKorendvke et alj|2006h 
aboard Hinode, high spectral and spatial resolution observa- 
tions in intense spectral lines formed at temperatures above 
1 MK have finally become available. In this work, by using 
combined SUMER (Solar Ultraviolet Measurement of Emitted 
Radiation)/SoHO and EIS/Hinode observations of the quiet Sun 
and a novel technique, we have characterised the Doppler shift 
versus temperature behaviour up to 2. 1 MK. 

In Section 2, we describe the SUMER data and their analysis. 
In addition, we determine the Doppler velocity for N v /11238, 
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Fig. 1. Fe XII 195 A EIT image. The areas observed on the disk 
by SUMER (smaller red box) and EIS (larger black box) are 
shown. They are located near disk centre in a quiet region. The 
off-limb pointing of the EIS slit (about 50 " above the west hmb) 
is indicated by the vertical black line. 



N V /11242, and Mg x A625 lines in the quiet Sun near disk centre. 
We also discuss the EIS/Hinode data, their analysis, and their 
co-alignment with the SUMER data. In Section 3 we present our 
method to measure the average Doppler velocity of hot coronal 
ions using SUMER/SoHO and EIS/Hinode data. Section 4 and 
5 present results and conclusions, respectively. 



2. Observations and data analysis 

Spectra of a quiet region near disk centre were taken on 6 April 
2007 from 00:00 to 06:00 UTC during the first EIS-SUMER 
campaign (HOP 3). Figure[T]shows the Sun on the day of obser- 
vations as recorded by the Extreme Ultraviolet Telescope (EIT, 
iDelaboudiniere et"al]|1995h aboard SOHO. The quiet Sun areas 
observed by SUMER and EIS are visible near disk centre. 



2.1. SUMER data: Doppler shift of Mg x and N v 

SUMER is a normal incidence spectrograph operating over the 
wavelength range 450 A to 1610 A. It is a powerful UV in- 
strument capable of making measurements of bulk motions in 
the chromosphere, TR, and low corona with a spatial resolu- 
tion of 1" across and 2" along the slit and a spectral scale 
of ^ 43 mA/pix at 1240 A (first order) (IWilhelm et al.ll 19951: 
iLemaire et al .119971) . which is high enough to measure line shifts 
down to 1 km s^'. Spectral images have been decompressed, 
wavelength-reversed, corrected for dead time, flat-field, and de- 
tector electronic distortion. 

Since SUMER does not have an on-board calibration source, 
we derive the wavelength scale by using a set of chromospheric 
lines assumed to have negligible or low average Doppler shift on 
the quiet Sun. In short, caUbrating takes three main steps: 

1. Identify the spectral lines by using the preliminary wave- 
length scale. 

2. Find the exact pixel position of the lines by fitting Gaussian 
curves. 

3. Perform a polynomial fit to find the dispersion relation. 
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To perform an accurate calibration we have to choose fairly 
strong and unblended lines from neutral and singly ionised 
atoms. For instance. Si i, S i, and C i transitions, occurring at 
temperatures between 6500 and 10000 K, are very suitable cal- 
ibration lines, since lines formed at su ch temperatur es show, on 
average, no substantial Doppler shifts jSamai 

The SUMER data analysed here consist of a raster formed by 
40 slit positions with a step size of 1.87" taken between 03:26 
and 04:31 UTC by exposing through the 1" x 120" sUt for 100 s. 
Recorded spectra cover the 1210 A to 1270 A range contain- 
ing the N v Unes at 1238.8 A and 1242.8 A (log(r/[K])=5.24) 
and the Mg x line at 625 A (log(r/[K])=6.00). The above spec- 
tral range also contains a forbidden line of Fe xii at 1242 A that 
would be ideal for a combined study with EIS. However, on the 
quiet Sun this line is w eak and blended with chr omospheric lines 
from Ca ii and S i (see lDammasch et"d]ll999ah and, for this rea- 
son, it is not possible to derive its average Doppler shift with the 
accuracy needed for our study. 

After checking that there are no significant instrumental 
drifts in the position of the spectral lines (either along the slit 
either in time across the raster) we have obtained a high signal to 
noise spectrum to be calibrated in wavelength. We have selected 
ten suitable reference lines in the part of the spectrum which con- 
tains the N V lines. A first-order polynomial fit of their positions 
in pixels versus their laboratory wavelengths, leaves residuals of 
less than ± 0.2 pixels (less than 2 km s"', see Fig.|2]). 

For the part of the spectrum containing the Mg x line, we 
have selected nine suitable reference lines. If we use a first-order 
polynomial fit to the line positions, the residuals up to ± 0.6 pix- 
els (top panel of Fig.[3]l. The residuals clearly show the need for 
a second-order polynomial fit that yields residuals of + 0.05 pix- 
els (bottom panel of Fig. O. The need for the higher order fit 
comes most likely from residuals in the correction of the elec- 
tronic distortion of the detector image. 

From the above analysis, for the N v /11238 and A\2A2 
lines we find net red shifts of (9.9 + 0.6) km s ' and (12.1 + 
0.6) km s"' by assuming the rest wavelengths to be 1238.800 A 
and 1242.778 A (Edlen 1934) . These velocities agree with those 
from lines formed at similar temperatures repor t ed in the liter- 
ature (e.g.,'Brekke et al. "T997I: IChaTet al.lll998t IPeter & Judgd 
[1999; Teriaca etal. 1999a). 

In the case of the Mg x A625 line (observed in the second- 
order of diffraction), we use 624.967 A as rest wavelength (aver- 
age between the values of (624.965 ± 0.003) A (Dammasc h et alJ 
[T999ah and (624.968 ± 0.007) A (iPeter & JudgellT999h '). obtain- 
ing a Doppler velocity of (-1.8 + 0.6) km s"^ (Table [1). With 
the above rest wa veleng th, t he Doppler s hifts measured in Mg x 
by iBrekke et all (119971) and IChae et alj (Il998h would result in 
velocities of about km s~ 1 and -2 km s~ 1, respectively. 

IBrekke et"an(fT997h discuss a possible blend in the red wing 
of the Mg X A625 Hne with a first-order Si 11 1250.089 A line and 
conclud that it has a negligible effect on the derived line shifts 
in spectra taken on t he KBr-coated part of the SUMER detector 
iTeriaca et al.l (l2002h also discuss the relevance of blending from 
first-order lines to the Mg x line in the quiet Sun, concluding that 
they account for about 10% to 15% of the signal in the observed 
profile when observed on KBr (mostly from the Si 11 1250.089 A 
line) and to 1% to 2% when observed on the uncoated detector 
(due to the low sensitivity to first-order lines with respect to the 
KBr-coated parts). Since our observations were obtained on the 
uncoated part of the SUMER detector, we consider blends from 
first-order lines as entirely negligible. 



Table 1. Doppler velocity of N v and Mg x lines: 



ions 


rest wavelength 


V 


6v 


(log (r/[K])) 


[A] 


[km s-i] 


[km s-'] 


N V (5.24) 


1238.800 


9.9 


0.6 


N V (5.24) 


1242.778 


12.1 


0.6 


Mg x (6.00) 


624.967 


-1.8 


0.6 



g 1000 r 

3 800 - 

m - 

g. 600 r 

"3 400 r 

■q. 200 \ 

1225 



N V calibration, first order 




1225 



1230 



1235 



1240 1245 
wavelength /A 



1250 



1255 



1260 



Fig. 2. Upper panel: first-order polynomial fit to pixel positions 
versus rest wavelengths of 10 suitable reference lines around the 
N V lines. Lower panel: fit residuals are under 0.2 pixels (less 
than 2 km s"'). 



Blends from second-order lines such as S x 624.70 A, O iv 
624.617 A, and 625.130 A are potentially more dangerous. 
However, these lines are far enough from the centre of the 
Mg X line and would produce detectable "shoulders" in the ob- 
served profile. Moreover, constrained double Gaussian fits show 
no significant improvement in line fitting over fits with a sin- 
gle Gaussian. This speaks against significant blending in the 
Mg X A625 line also from second-order lines. 

2.2. EIS data: analysis and co-alignment 

The EIS instrument on Hinode produces high-resolution stig- 
matic spectra in the wavelength ranges of 170 to 210 A and 250 

to 290 A. The instrument has 1 spatial pixels and 0.0223 A 
spe ctral pixels. More details are given by ICulhane et alj (l2007h 
and lKorendvke et al.1 (|2006). 

Our EIS data consist of four raster scans and one sit-and- 
stare sequence. The four scans were taken at disk centre on the 
same target area as observed by SUMER, whereas the sit-and- 
stare sequence was obtained above the limb. EIS rasters are ob- 
taining by scanning from west to east, in the opposite direction to 
the SUMER data discussed here. Raster 4 is the closest in time 
to the SUMER scan. Data were reduced by using the software 
provided within SolarSoft. 

After accounting for the offset between the short and the long 
wavelength ranges of EIS, we use pairs of radiance maps ob- 
tained in lines from similar (at least not too different) tempera- 
tures to perform a co-alignment between scans made by the two 

instruments with an accuracy of roughly 1 (see Fig.|4li. The em- 
ployed line pairs include N v ^1238 (log(r/[K])=5.24, SUMER) 
and He 11 A256 (log(r/[K])=4.7, EIS), as well as Mg x ^625 
(log(r/[K])=6.00, SUMER) and Fe x ^184 (log(r/[K])=6.00, 
EIS). 

Since the wavelength range covered by EIS does not contain 
any suitable neutral or singly ionised line (He 11 A256, besides be- 
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Fig. 3. Top pair of panels: first-order polynomial fit to pixel posi- 
tions versus rest wavelengths of 9 suitable reference lines around 
the Mg X lines. Residuals are large, up to 0.6 pixels. Bottom pair 
of panels: second-order polynomial fit to pixel positions versus 
rest wavelengths. Residuals lie within + 0.05 pixels. 



Fig. 4. Radiance maps of the common area of study in 
N V ^1238 A (SUMER, left) and He n A256 A (EIS, right). 
Contours of He ii A256 are plotted on both maps. The red boxes 
show the area in common to SUMER and the short and long 
wavelength ranges of EIS. 



and Si x that also have diff'erent mass-to-charge ratios (5.1 and 
3.1), although the difference is smaller. 

It should also be noted that Fe and Mg have very similar 
first ionisation potentials (FIP), excluding eventual differences 
in the average flows due to the fraction ation pro cesses believed 
to occur in the chromosphere (see, e.g.. lLaming|i2004) . 

From the above considerations, we conclude that assuming 
the average velocity of Fe x ions to be equal to that of Mg x ions 
is justified and should not result in a systematic error that would 
be large enough to change our results and conclusions substan- 
tially. Therefore, in the following, we take Fe x /1184 as our ref- 
erence line with respect to which the average Doppler shift of 
the remaining spectral lines recorded by EIS is determined. 



ing optically thick, is formed at too high a temperature to be un- 
shifted), an absolute wavelength calibration is not possible using 
the EIS data alone. We, therefore, take the approach of finding 
spectral lines recorded by the two instruments that are formed 
roughly at the same temperature by assuming that their average 
quiet Sun velocities are the same. Specifically, we assume that 
the average velocities of Mg x and Fe x ions are the same in the 
common observed area and equal to (-1.8 ± 0.6) km s"'. 

The above assumption requires some further consideration 
as it cannot be immediately excluded that the different ions have 
different velocities. In fact, mechanisms that can selectively heat 
and accelerate ions, such as ion-cyclotron resonance absorption 
of Al fven waves (see, e.g.. iTu & Marschlll997l: iMarsch & ful 
Il997h . have long been considered as candidates for the solar 
wind acceleration. This mechanism is usually considered for 
low-density, open fiel d regions such as coronal holes. However, 
iPeter&Vocksl (120031) claim that it may also operate at the base 
of open funnels in the quiet Sun, some 2 to 10 Mm above the 
photosphere. Since Fe x and Mg x have substantially different 
mass-to-charge ratios (6.2 and 2.7, respectively), the Fe x ions 
could be accelerated more (due to the lower gyro-frequency) and 
at lower heights (due to the expected fall-off of the magnetic field 
strength with height) than the Mg x ions. Thus, we may be un- 
derestimating the Fe x average outflow speed by taking it to be 
equal to that of Mg x. On the other hand, at slightly higher tem- 
peratures, we obtain very consistent results from lines of Fe xii 



3. A method of measuring the Doppler velocity of 
EIS coronal lines 

Above the limb we assume that spectral lines have, on aver- 
age, no Doppler shift (motions out of the plane of sky cancel 
out on average). This is verified by plotting histograms of veloc- 
ity (with respect to the average) and histograms of the distance 
in wavelength between a line and the reference Fe x /1184 line. 
These histograms are very narrow, nearly Gaussian, and have an 
FWHM of 3 km s"' to 4 km s"' (Figs. |5] IS. 

From the above assumption^Q, the average of the off-limb 
line-distance histogram is then 



Ref 



(1) 



To be precise, the quantity AAoff we measure from the histogram of 
the distance in (off-limb) wavelength should be corrected for the effect 
of solar rotation before being used at disk centre. This requires the right- 
hand side of Eq.[T]to be multiplied by 1 + v^^Jc. In the case of the Sun, 
I'lot ~ 2 km s"' and the term is practically equal to one. 

- Also the gravitational red shift plays no role here. In fact, with 
EIS we are comparing distances between spectral lines on the disk and 
above the limb, and the effect, because the gravitational red shift is the 
same at the two places, cancels out. In the case of the SUMER data, the 
Doppler shift are measured with respect to the chromosphere, which is 
assumed to be at rest. 
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where A^"^^ and are the rest wavelengths of the reference line 
and of the line whose shift we want to measure. 

On the disk we measure the quantity l^A - A - A'^'^f at each 
pixel or group of pixels (where /l*''^ is the wavelength of the 
reference line). This can be written as 



AA^ A0 + 6A - aI'^ - SA'^'f = ^Aoff +6A - 6A'^'^. 



(2) 



where 6 A and dA'^"^ are the line Doppler shifts. Considering that 



V = c^, we have 



AA = AAoff + - ^0 - 



(3) 



If we call 6v = V - v**"^, the difference between the velocities 
measured in the two lines, then 6v becomes 



dv = — \AA - AAoff 1 1 + 

Ao\ -"1 c 



(4) 



Considering now the average value of v'^"^, v**"^, we can write 



AA - AAoff 1 1 + 



(5) 



We measure AA by taking the average of the distribution of 
wavelength differences between the two lines . Then, using equa- 
tion (|5]l we calculate the relative velocity dv of hot coronal ions 
with respect to Fe x /1184 that we assume to have an average 
Doppler velocity v*''^ equal to the average Doppler velocity of 
Mg X A625 (-1.8 + 0.6 km s"'). 

Because of thermal variation along the 98.5 min Hinode or- 
bit, there is a shift in the spectral line positions with time that 
needs to be corrected before using EIS data. Routines for cor- 
recting this variation are provided with the data reduction pack- 
age and the correction has been applied to our data. However, we 
are not sure that the above correction (based on the data itself) 
is reliable down to the accuracy needed for the measurements 
discussed here («; 2 km s"'). A major advantage of the tech- 
nique presented in this paper is that we do not rely on this cor- 
rection because we are just measuring the spectral distances be- 
tween spectral lines0. We only need to make sure that all spectral 
lines move across the detectors in the same way, i.e., that there 
are no changes in the image scale with time. The top panel of 
Fig. |7] shows the dependence on time (and solar X) of line posi- 
tions (averages along the slit, in wavelength units) for three lines 
(Fe X /1 184, Fe x A257, and Fe xii /1 193) when no correction for 
temporal drifts is introduced. Offsets in wavelength (indicated in 
the figure) have been subtracted to make the comparison more 
evident. The behaviour of all the lines is very similar. The bot- 
tom panel of the figure shows line distances (differences in the 
values in the top panel) versus time and solar X. It can be seen 
that changes in line distance (wavelength differences) are very 
small (sil mA, « 1.5 km s"') and show no significant depen- 
dence with time. Hence our assumption is fulfilled. 

To improve the signal-to-noise ratio, before fitting the line 
profiles we have binned the EIS on-disk spectra over three raster 
positions and over three pixels along the slit (3x5 binning for the 
Fe XV /J284, weak on the quiet Sun). A 2x2 binning was used for 
the off-limb spectra. 
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Fig. 5. The histograms of line position in wavelength for 
Fe XII /1 193 line off-limb (dotted line) and at disk centre (solid 
line). The off-limb histogram is narrower than the on-disk one 
because net mass motions are essentially absent. 
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Fig. 6. The histogram of distance in wavelength between 
Fe XII /1 195 line and the reference line (Fe x /1 184) off-limb. The 
histogram is nearly Gaussian and with a sigma of ^ 1.3 km s"', 
that can be taken as an estimate of the instrument accuracy in 
measuring the distance in wavelength. 



4. Results 

Using the technique discussed in Section|3]we have obtained the 
line of sight velocities from the Doppler shifts of many coro- 
nal lines formed at temperatures between 1 MK and 2 MK in 
the quiet Sun. These lines, their deduced relative (to Fe x /1 184) 
Doppler shifts and the corresponding uncertainties are listed in 
Table |2l In this table, we report the values obtained in the four 
EIS raster scans together with the estimated errors. Errors are ob- 
tained by error propagation analysis of Eq.|5] where the error on 
AA^)f f is given by the width cr of the off-limb line-distance distri- 
bution. On the other hand, the on-disk line distance distribution 
is also broadened by the different Doppler speeds characterising 
the two lines at different locations. For this reason the error for 
AA is assumed to be equal to that of AAoff, leading to the same 



^ This method also ensures that there is no bias from the observed 
line intensities in the derived average velocity. 



5 



Dadashi, Teriaca, Solanki: The quiet Sun average Doppler shift of coronal lines 



0.04 



S 0.02 



0.01 



0.00 



solar X (arcsec) 
10 20 

iis Fe XIII M93 (\-193.50) 
O Fe X XS67 (X-257.27) 
A Fe X \184 (\- 184.54) 



* M * X 



> 

A _ 



A A 



i AAA' 



3000 2500 2000 1500 1000 
Time (sec) 

solar X (arcsec) 

-10 10 20 



500 



A Fe XIII M93 - Fe x M84 
o Fe XIII X193 - Fe X 



A A A 



i " A i 

^^a'^ AA^A^ 



0.012 - 



3000 2500 2000 1500 1000 
Time (sec) 



Fig. 7. Upper panel: line positions, in wavelength units, of the 
Fe XII /1193, Fe x /1257, and Fe x /1184 lines versus time (and 
solar X) at disk centre. The wavelength values indicated in the 
panel have been subtracted to ease comparison. Lower panel: 
wavelength difference between spectral lines for the data in the 
upper panel versus time (and solar X). It can be seen that changes 
in line distances (wavelength differences) are very small (wl 
mA, 1.5 km s ') and show no dependence on time. 



error for a given line in all the rasters. An independent estimate 
of these errors can be obtained by looking at the standard devi- 
ation of the measured values of 6v obtained in the four rasters. 
From Table|2] values of 0.1 km s"' and 0.7 km s"' would be ob- 
tained for Fe XII A195 and Fe xv /1284, respectively. These values 
are lower than our estimated errors and indicate the most likely 
we do not underestimate the errors provided in Table [2] 

In the following, we take the average of rasters # 3 and 
# 4 as representative values of the average Doppler shift on the 
quiet Sun, the closest in time to the SUMER observations. By 
adding the value of average Doppler velocity of our reference 
line (Mg x A625) to these values (v - v^'^^ + 6v), the absolute 
Doppler velocity of lines are calculated (see TableO. The error 
reported in this table, 5(y), is obtained by summing in quadrature 
the errors on relative line Doppler shifts (reported in Table |2]i, 
and the error in measuring the average Doppler velocity of our 
reference line (Mg x A625). The values in Table O are plotted in 
Fig.[8]and as a blow-up for temperatures above 1 MK in Fig.|9] 

Figure |8] displays the Doppler shift in the quiet Sun at disk 
centre of various TR and coronal ions in SUMER and EIS spec- 
tra. Positive values indicate red shifts (downflows), while nega- 
tive values indicate blue shifts (upflows). Values obtained in this 
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Fig. 9. Magnification of the high-temperature part of Fig. [8] 
showing the results obtained in this work for coronal lines. 
Symbols and lines as in Fig. [8] 



investigation (both SUMER and EIS) are marked by stars. In the 
figure the results of , Teriaca et al.. (199 9a^), Peter & Judge, (19991) . 
and the theoretical Doppler shifts of iHansteen et al.l ( 20101) are 
also shown. 

The solid and dashed lines respectively represent a by-eye 
fourth and sixth order polynomial representation of the com- 
bination of our pres ent work for TR and coronal lines (T < 
2 MK) and those of iTeriaca et alj (Il999ah . Thus, the v versus 
log(r/[K]) behaviour can be parametrised by the coefficients 
that are listed in Table |4] 

At temperatures higher than 1 MK our results show a mod- 
est but clear increase of the average blue shift with respect to the 
value measured around 1 MK (-1.8 km s"'), reaching a maxi- 
mum of about -4.4 km s"' at the Fe xiv formation temperature 
of 1.8 MK. At even higher temperatures, the average line shift 
appears to decrease again, reaching a value of (-1 .3 +2.6)km s"' 
obtained in Fe xv A284 line at 2.1 MK. 

To check if our results are dependent on the choice of our 
reference line (Fe x /il84), we have used another Fe x line 
(/i257) and calculated all Doppler velocities with respect to it. 
The Doppler velocities are found to agree within 1 km s ' with 
those in Table |2] 

It is clear that the absolute values of the velocities derived 
here depend on the value adopted for the reference line (in this 
case -1.8 km s '). To check the dependence of our results on 
the value of the velocity of our reference line (Fe x /II 84, which 
we have supposed to have the same average Doppler velocity 
of Mg X A625), we used Eq. |5]to calculate the absolute veloc- 
ity for different values of the rest wavelength of the Mg x line 
(which ultimately determine the velocity in the reference line). 
The dependence of the Doppler shift on the Mg x rest wave- 
length and on the resulting Doppler velocity is shown in Fig.fTOl 
for the Si x ^261, Fe xii .1195, Fe xiii .1251, Fe xiv A264, and 
Fe XV /1284 l ines. It sho uld be noted that using the rest wave- 
length from iKellvl (11987 ). would result in red shifts for most of 
the coronal lines studied here. 



5. Summary and conclusions 

In the present study we measured the Doppler shift of TR and 
coronal lines on the quiet Sun (see Tables 1 & 3). We measured 
the absolute average velocity of Mg x A625, N v /11238, and 
N v /1 1242 lines using SUMER data obtaining a blue shifted 
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Fig. 8. Average Doppler shift in the quiet Sun at disk centre of various TR and coronal ions measured from SUMER and EIS 
spectra. Positive values indicate red shifts (downflows) while negative values indicate blue shifts (upflows). Solid and dashed lines 
respectively represent fourth and and sixth order polynomial curves by-eye fitted to the combination of the measurements obtained 
in this work and those of .Teriaca et al., (,1999a) . 



Doppler velocity of (-1.8 + 0.6) km s ' for Mg x A625 and red 
shifts (downward motions) of (9.9 + 0.6) and (12.1 + 0.6) for 
the N V /11238 and A\2A2 lines. Our result s for TR lines agree 
with pre vious stu dies done with SUMER (Brekke et al."l997^, 
IChae et al. 1998; Peter & Judgg 1999: Teriaca et aL, 1999a) . Our 
result for the Mg x line is also consistent with the literature, once 
the differences on the adopted rest wavelength are accounted 
for (see Section 2.1). Of course, our value for the average shift 
in the Mg x line depends on the adopted rest wav elength of 
624.967 A ( IDammasch etal]|1999al: [Peter & Judgel lT9"99. av- 
erage between the values of). Since the velocities obtained from 
the EIS lines also depend on the absolute velocity of the Mg x 
ions, the precise determination of the rest wavelength of Mg x 
is paramount. The accurate, indepe n dent, a nd hi ghly consistent 
measurements of IDammasch et al] (Il999ah and 'Peter & Judg j 
119991) are the best available so far. In any case. Fig. [TOl allows 
assessing the effect of a different rest wavelength for the Mg x 
A625 Une on the resulting Doppler shifts. 

By using cospatial and nearly cotemporal SUMER & EIS 
data we, for the first time, measured the average Doppler shift of 
several coronal lines with formation temperatures above 1 MK 
in the quiet Sun near disk centre. We obtained negative Doppler 
shifts (upward motions) between -2 to -4 km s"' at temperatures 
between 1 and 1.8 MK. At the highest temperature we can in- 
vestigate, T - 2.1MK (Fe xv /1284), we find indications of a 
reduction of the observed blue shift down to 1.3 km s"', a value 



that is compatible with no line shift, but also with a blue shift of 
up to 4 km s"^ 

As discussed in the Introduction, 1-D models assuming en- 
ergy release in magnetic loops had successfully managed to re- 
produce the observer red shif t at TR temper atures and blu e shifts 
at co ronal temperatures (e.g. jTeriaca et al.f l999b; Soadar oet al.l 
2006). However, these models do not explicitly consider the 
heating mechanism. 

More recently, comprehensive 3-D models of the solar at- 
mosphere have been developed. Models assuming that coronal 
heating is caused by Joule dissipation of currents produced by 
stressing and braidi ng of the rnagnet i c fiel d produce red shifts 
at all temperatures dPeter et alj|2004 2006). which is con trary 
to the present result s and those of dPeter & Judge! 1 199^ and 
dTeriaca et alJll999ah . 

The addition to the model of episodic injections of emerg- 
ing magnetic flux that reconnects with the existing field, pro- 
duces rapid, episodic heating of the upper chromospheric plasma 
to coronal temperatures (Hansteen et al. 2010). Doppler shifts 
are computed by these authors for a model with an unsigned 
magnetic field of | B, |^ 75 G, which is compatible with the 
quiet Sun (Trujillo Bue no et alj2004tlDonu'nguez Cerdena et al.l 
2004 iDanilovic et al.l l2010l) . The resulting Doppler shifts are 
roughly consistent with those reported here at all temperatures 
up to the Fe xiv formation temperature (log(r/[K])=6.26) (see 
Fig. [8]). Above this value, the line shifts predicted by the model 
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Table 2. Doppler velocity of ions with respect to Fe x /1184. 



line 6v 6(6v) 







[km s 






[km s"'] 




# 1 


#2 


#3 


#4 




Fe X 257.262 


0.86 


0.28 


0.52 


1.44 


2.30 


Fe XI lsU.4Ul 


0.96 


1 1/1 

1.14 


0.76 


1.22 


2.11 


Fe XI loo. 216 


0.15 


ATI 

-0.71 


A C 

-0.25 


-0.68 


1 oo 

1.88 


Fe XI 188.299 


-1.35 


-2.35 


-1.41 


-1.57 


2.14 


Si X 258.375 


-0.41 


A OT 

-0.87 


-1.57 


1 1 A 

-1.10 


2.04 


Si X 261.058 


-0.96 


-0.04 


-1.19 


-0.73 


2.29 


Six 271.990 


-1.51 


-1.62 


-2.06 


-2.17 


2.07 


S X 264.233 


0.21 


-0.01 


-0.47 


-0.47 


2.09 


Fe XII 193.509 


-0.64 


-0.83 


-0.16 


-0.46 


1.98 


Fexii 195.119 


-0.52 


-0.67 


-0.52 


-0.67 


1.52 


Fe XII 256.925 


-0.54 


-1.00 


-1.70 


-0.65 


2.76 


Fe XIII 202.044 


-0.43 


-0.73 


-0.14 


-0.43 


2.02 


Fe XIII 251.953 


-1.54 


-3.08 


-2.84 


-1.42 


2.96 


Fe XIV 264.787 


-2.62 


-2.96 


-2.27 


-2.84 


2.13 


Fe XIV 274.203 


-2.47 


-2.14 


-3.45 


-1.59 


2.13 


Fexv 284.160 


1.38 


0.33 


1.28 


-0.31 


2.56 



Columns 2 to 5 list the results for the four different rasters. 
The last column provides the error, which is equal for all 
rasters (see text). 




0.9B 0.97 0.96 0.95 

^R.f Mg X ( + 624 ) 

I , , , , I , , , , I , , , , I , , , , I 

-10 -5 5 10 

^Ref ^ ^'^^ ^ {]im s"') 



Fig. 10. Variation in the absolute Doppler velocity of hot coro- 
nal ions as a function of the rest wavelength of the Mg x line and, 
consequently, of the Doppler velocity from the Fe x /1184 refer- 
ence line. The vertical dotted line corresponds to the adopted rest 
wavelength of 624.967 A. 

return to zero and become slightly red shifted at the Fe xv for- 
mation temperature (log(r/[K])=6.32). 

Our Fe xv observations, instead, yield a blue shifted value. 
However, the error bar is large enough that no firm conclusion 
can be drawn. A real discrepancy at this temperature between 
the observations and the model may possibly be due to bound- 
ary effects in the adopted computing box or, more simply, to the 



Table 3. Absolute average Doppler velocity of ions (averaged 
over rasters #3 and #4). 



line uogW IV'^i)) 


1 T rii^m c ^1 
V [Jvlll a J 


scv^ rvm Q^h 


Fex 184.536 (6.00) 


-1.8'' 


0.6'' 


Fe X 257.262 (6.00) 


-0.82 


2.38 


Fexi 180.401 (6.04) 


-0.81 


2.83 


Fexi 188.216(6.04) 


-2.27 


1.97 


Fexi 188.299 (6.04) 


-3.29 


2.22 


Six 258.375 (6.15) 


-3.14 


2.13 


Six 261.058 (6.15) 


-2.76 


2.37 


O' '^Tl AAA / IZ 1 C\ 

Si X 271.990 (6.15) 


-3.92 


2.16 


S X 264.233 (6.15) 


-2.27 


2.17 


Fe XII 193.509 (6.15) 


-2.11 


2.07 


Fexii 195.119(6.15) 


-2.40 


1.63 


re XII zjo.yzj (0. 1 j) 


-/.vo 




Fe XIII 202.044 (6.20) 


-2.09 


2.11 


Fe XIII 251.953 (6.20) 


-3.93 


3.02 


Fe XIV 264.787 (6.26) 


-4.36 


2.21 


Fe XIV 274.203 (6.26) 


-4.32 


2.21 


Fexv 284.160 (6.32) 


-1.34 


2.63 



^ The last column provides the error on the absolute average 
Doppler velocities (see text). 
Absolute velocity and error assumed equal to that of Mg x. 



Table 4. Coefficients of the fourth and sixth order polynomial 
representations of the average Doppler velocity versus tempera- 
ture behaviour of the quiet Sun. 



coefficient 


4'* order 


6''' order 


ao 


7874.8 


-185202.25 


ai 


-6427.5 


223219.47 


a2 


1943.26 


-111257.36 


as 


-257.766 


29345.005 


a4 


12.6596 


-4319.1993 


as 




336.35535 


ae 




-10.827987 



All decimals shown here are significant. 

fact that the turning point to zero or red shifted values occurs in 
reality at temperatures even higher than those measurable in this 
study. However, it may also be possible that important ingredi- 
ents are still missing in the theoretical/numerical description. 

Hans teen et al.l (1201 Oh suggest that the high speed flows 
generated by the localised heating events may be related to 
the high speed upflows that have been deduced from sig- 
nificant blueward asymmetr ies in Hinode/EIS observations of 
many TR and coronal lines (Hara et al."2008t iDe Pontieu et al.l 
[2009; Mcintosh & De Pontieu 2009b a) that have been linked to 
spicules ((De Pontieu et al.ll201 II) . Thus, a more accurate intro- 
duction of the effect of spicules into the model may be necessary. 

The model also foresees a smaller amount of Doppler shifts 
at both TR (lower red shifts) and coronal temperatures (lower 
blue shifts). However, this is not surprising (as the model does 
not specifically attempt to reproduce the solar region studied 
here), and the model is most likely compatible with our mea- 
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surement. Since the model does not consider open field lines, 
this would indicate that the observed line shifts are simply a 
by-product of the heating of, and related mass supply to, closed 
loops. This would be consistent with the finding that the plasma 
outflows seen on the quiet Sun, such as the Ne vm blu e shift, are 
largely associated with closed loops dTian et al.ll2008h . 

However, it cannot also be excluded that hot open funnels 
significantly contribute to the blue shifted emission observed be- 
tween 1 and 2 MK. Only a detailed analysis of spatially resolved 
flows compared to field extrapolations (similar to what was done 
for the Ne vm line at lower temperatures by Tian et af] l2008h 
would shed some light on the question. Given that these hotter 
lines are weak on the quiet Sun, this analysis is very difficult and 
will require a VUV spectrometer with much higher throughput 
than currently available, such as the spectrometer proposed for 
the Japanese Solar C mission (Teriaca et al. 201 1). 
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